Identification of a polymer growth process with an equilibrium multi-critical collapse 
phase transition: the meeting point of swollen, collapsed and crystalline polymers 
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We have investigated a polymer growth process on the triangular lattice where the configurations 
produced are self-avoiding trails. We show that the scaling behaviour of this process is similar to 
the analogous process on the square lattice. However, while the square lattice process maps to the 
collapse transition of the canonical interacting self-avoiding trail model (ISAT) on that lattice, the 
process on the triangular lattice model does not map to the canonical equilibrium model. On the 
other hand, we show that the collapse transition of the canonical ISAT model on the triangular lattice 
behaves in a way reminiscent of the S-point of the interacting self-avoiding walk model (IS AW), which 
is the standard model of polymer collapse. This implies an unusual lattice dependency of the ISAT 
collapse transition in two dimensions. 

By studying an extended ISAT model, we demonstrate that the growth process maps to a multi- 
critical point in a larger parameter space. In this extended parameter space the collapse phase 
transition may be either S-point-like (second-order) or first-order, and these two are separated by a 
multi-critical point. It is this multi-critical point to which the growth process maps. Furthermore, 
we provide evidence that in addition to the high-temperature gas-like swollen polymer phase (coil) 
and the low-temperature liquid drop-like collapse phase (globule) there is also a maximally dense 
crystal-like phase (crystal) at low temperatures dependent on the parameter values. The multi- 
critical point is the meeting point of these three phases. Our hypothesised phase diagram resolves 
the mystery of the seemingly differing behaviours of the ISAW and ISAT models in two dimensions 
as well as the behaviour of the trail growth process. 
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I. INTRODUCTION 

Over the past 25 years various lattice models of a single 
self-interacting polymer chain in (dilute) solution have 
been analysed in both two and three dimensions. The 
fundamental physical phase transition [1] that these mod- 
els attempt to mimic is that of the collapse of a single 
polymer in a poor solvent as the temperature is lowered. 
At high temperatures a polymer is swollen relative to a 
reference Gaussian state, while at low temperatures the 
polymer forms a liquid-drop-like globule [r,'2;. There is a 
continuous phase transition expected between these two 
states, which is referred to as the 6'-point. One question 
that arises concerns the robustness of the universality 
class of the collapse transition. The standard theory f3|- 
[S] of the collapse transition is based on the n — > limit 
of the magnetic tri-critical 0"* — cf)^ 0{n) field theory and 
related Edwards model with two and three body forces 
[51 [7] , which predicts an upper critical dimension of three 
with subtle scaling behaviour in that dimension. On the 
other hand recent studies [SHIO] of semi-flexible polymers 
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indicate that collapse transition can be first-order, and 
a polymer crystal can exist at low temperatures in the 
presence of stiffness. 

The canonical lattice model of the configurations of 
a polymer in solution has been the self-avoiding walk 
(SAW) where a random walk on a lattice is not allowed 
to visit a lattice site more than once. Self-avoiding walks 
display the so-called excluded volume effect where they 
are swollen in size relative to unrestricted random walks 
of the same length: their size scales with a different char- 
acteristic exponent to that of the unrestricted random 
walk. A common way ^ to model intra-polymer inter- 
actions in such a walk is to assign an energy to each 
non-consecutive pair of monomers lying on neighbouring 
lattice sites. This is the interacting self-avoiding walk 
(ISAW) model, which is the standard lattice model of 
polymer collapse using self-avoiding walks. 

On the other hand self-avoiding trails [TT|, which are 
lattice random walks that are not allowed to visit a lat- 
tice bond more than once also display the same excluded 
volume behaviour as SAW and physical polymers, with 
the same scaling exponents. Moreover, a self-interacting 
self-avoiding trail model (ISAT), where multiply occu- 
pied sites are assigned an energy, has also displayed some 
of the characteristics of the polymer collapse described 
above P^15j . However, analyses of both two- and three- 
dimensional self-interacting trails p"51 - [T5] indicate that 
the collapse transition of the ISAT model is in a different 
universality class to that of self-interacting self-avoiding 
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walks (ISAW) in those respective dimensions. There is 
no clear understanding of why this is the case if true. 
It is important to note that work in two dimensions has 
focused on the square lattice model only. 

Some of the work [T31 [2] on the IS AT model has been 
via a growth process known as "kinetic growth trails" or 
"smart kinetic trails" that maps to one particular tem- 
perature of the equilibrium ISAT model on the lattices 
studied (it is important to note that these lattices are 
of coordination number four). It was proposed in |16j 
that the collapse transition associated with "smart ki- 
netic trails" is first-order. Clear evidence was produced 
in [16| to demonstrate that there was a first-order transi- 
tion in three dimensions on the diamond lattice. On the 
other hand no evidence of this could be found in two di- 
mensions. In fact, on the square lattice recent work |15| 
verified that the ISAT collapse transition on the square 
lattice is precisely that given by the growth process and 
is not first-order or like the 0-point of ISAW. 

Another approach has been to analyse a generalised 
ensemble (rather than the finite length ensemble) of the 
model via transfer matrix jl7j : this has produced some 
intriguing results with different values of critical expo- 
nents being estimated. It was pointed out in |17) that 
these results are not incompatible with those for the fi- 
nite length ensemble. This may be an indication though 
that the point in question is multi-critical of some type. 

In this paper we study the ISAT model and kinetic 
growth trails on the triangular lattice, importantly a co- 
ordination number six lattice. We demonstrate that the 
kinetic growth model does not map to any temperature 
of the canonical ISAT model but rather to a point (we 
shall call this the kinetic growth point) in the param- 
eter space of a generalised model we call the Extended 
ISAT model (elSAT). As such we have studied this elSAT 
model and identify this kinetic growth point as a multi- 
critical collapse point. By studying the elSAT model we 
build a picture of the collapse transition in the trail model 
of polymers that includes both the first-order transition 
suggested in |TS] and the ISAW f?-point. It may also 
prove useful in explaining the transfer matrix results. 

A. Review of previous results 

The collapse transition can be characterised via a 
change in the scaling of the size of the polymer with tem- 
perature. It is expected that some measure of the size, 
such as the radius of gyration or the mean squared dis- 
tance of a monomer from the end points, i?^(T), scales 
at fixed temperature as 

Rl{T)^An''^ (1.1) 

with some exponent v. At high temperatures the poly- 
mer is swollen and in two dimensions it is accepted that 
J/ = 3/4 jlSj. At low temperatures the polymer becomes 
dense in space, though not space filling, and the exponent 
is 1/ = 1/2. The collapse phase transition is expected to 



take place at some temperature Tj. If the transition is 
second-order, the scaling at Tt of the size is intermedi- 
ate between the high and low temperature forms. In the 
thermodynamic limit there is expected to be a singular- 
ity in the free energy, which can be seen in its second 
derivative (the specific heat). Denoting the (intensive) 
finite length specific heat per monomer by c„(T), the 
thermodynamic limit is given by the long length limit as 

C(T) = lim c„(T) . (1.2) 

One expects that the singular part of the specific heat 
behaves as 

C{T) ^ B\Tt~T\-" , (1.3) 

where a < 1 for a second-order phase transition. The 
singular part of the thermodynamic limit internal energy 
behaves as 

U{T) ^ B\Tt-T\^-°' , (1.4) 

if the transition is second-order, and there is a jump in 
the internal energy if the transition is first-order (an ef- 
fective value of a = 1). 

Moreover one expects crossover scaling forms [T^ to 
apply around this temperature, so that 

Cn{T)^n'"f'C{{T-Tt)n^) (1.5) 

with < < 1 if the transition is second-order and 

Cr^iT) ^ n CUT - Tt)n) (1.6) 

if the transition is first-order. From |19j we point out 
that the exponents a and cf) are related via 

2-a=i. (1.7) 



a. -point ISAW collapse The work of Duplantier 
and Saleur [20] predicted the standard 0-point behaviour 
in two dimensions, which has been subsequently verified 
[2T] . It is expected that 

= 3/7^0.43 and a = -1/3 . (1.8) 

Note that this implies that the specific heat does not 
diverge at the transition. However, the third derivative of 
the free energy with respect to temperature will diverge 
with exponent (1 + a)(j) ~ 2/7. At T = Tt it is expected 
that 

Rl{T) - ^n^/y g) 

b. ISAT collapse on the square lattice Previous 
work [15] on the square lattice has shown that there is 
a collapse transition with a strongly divergent specific 
heat, and the exponents have been estimated as 

= 0.84(3) and a = 0.81(3) . (1.10) 

At T = Tt it was predicted [H] that 

Rl{T)^An{\ognf . (1.11) 
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II. TRAIL GROWTH ON THE TRIANGULAR 
LATTICE 

Consider a stochastic process defined on the triangu- 
lar lattice as follows: Starting at an origin site, a lattice 
path is built up step-by-step by choosing between avail- 
able continuing steps from unoccupied lattice bonds with 
equal probability. The configuration produced is a self- 
avoiding trail or trail for short, where sites may be visited 
multiple times but bonds of the lattice are either visited 
once or not at all. If we do not consider the original oc- 
cupation of the origin as a visit, then each site may be 
visited up to three times on the triangular lattice. Note 
that when the process revisits the origin, the growth rule 
is slightly altered: the process may choose the direction 
of the first step as one of its options equally with the un- 
occupied lattice bonds. If it does this, a loop is formed 
and the growth process terminates. Apart from the ini- 
tial step, where there are six available steps, the number 
of available steps is therefore five minus twice the num- 
ber of previous visits. A picture of these three situations 
is provided in Figure [T] along with the probabilities of 
the next step that is added. In Figure [2] an example 
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FIG. 1; This figure illustrates the growth process on the tri- 
angular lattice. Steps are created with difi'erent probabilities: 
the probabilities of the next steps (not shown) are displayed 
underneath each case. Three general possibilities exist Eis the 
trail grows: firstly, that the trails enters a site that is previ- 
ously unoccupied — in that case the five choices of next step 
are chosen at random with equal probability, that is, 1/5; 
secondly, that the trail enters a site that has been previously 
visited once — in that case the three remaining choices of 
next step are chosen at random with equal probability, that 
is, 1/3; finally, that the trail enters a site that has been previ- 
ously visited twice — there is only one choice of continuation 
and this is made with probability 1. 

of a configuration produced by the process is illustrated. 
Apart from the origin, each singly visited site contributes 
a factor (g) to the overall probability, while each twice 
visited site contributes a factor (|)(|) and each triply- 
visited site contributes a factor (i)(i)(l). This model is 
the triangular lattice version of the growth model con- 
sidered previously [TBI. 

We will not count the initial occupation of the origin 
as a visit. If we denote the number of steps of the trail as 
n, and the numbers of singly, doubly and triply-visited 
site as mi, m2 and 7713, respectively, then these satisfy 




FIG. 2: An example of a trail with 13 steps on the tri- 
angular lattice. This trail has six singly visited sites, 
two doubly-visited sites and one triply-visited site. Note 
that the path may cross or touch at a visited site. This 
trail is produced by the growth process with probability 

- the individual 

fractional probabilities are the step probabilities in order of 
creation of the steps. 

The probability of a configuration, (p„, of n steps is 
denoted pci^n), and we have for each n that 

Y^PciV^n) ^ Pn , (2.2) 

where Pn is the probability that the growth process 
reaches length n. Let us define the expectation values 
of the number of singly, doubly and triply-visited sites 
per unit length ej{n) for the growth process conditioned 
on the process making it to length n respectively as 

ej{n) = = —^"^PGiVn) mj{ipn) , (2.3) 

and the fluctuations in these as fj{n) with 



Note that from Equation ( |2.1[ ) we have 

ei(n) + 262 (n) + 363(71) = 1 . (2.5) 

In the square lattice case [351 123] this stochastic model 
can be mapped to a specific temperature of an equilib- 
rium model, and moreover [131 115) it was shown that 
this temperature was a critical point of the equilibrium 
model. There, the exponents a and <j) defined in the in- 
troduction can be related to the behaviour of 62 and /2 . 
Assuming a similar form of critical behaviour is found in 
the triangular lattice model we would expect that 

^2(")-^^2-^^(3^ (2.6) 

and 

e3(-) - ^3 - ■ (2.7) 

Also 



n = mi + 2m2 -|- 3m3 



(2.1) 



(2.8) 
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and 



+ F3 



(2.9) 



f3(n) 



where either the n°"^ or the constant term dominates de- 
pendent on whether a is positive or negative. 

We have simulated this growth process for lengths up 
to n = 2^0 = 1,048,576 producing 7,120,000 samples 
calculating estimates for 62(71), 63(71), f2{n), and fz{n). 
In Figure [3] we plot f2(n) and 73(71) against 71: we im- 
mediately note that both quantities diverge. Taking the 
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FIG. 3: Plot of the fluctuations in the number of doubly and 
triply- visited sites f2{ji) and /3(n) against n. The slope gives 
us an estimate of acj). 

better behaved data for /^{n), we have estimated 

a0 = 0.734(6) . (2.10) 
We can therefore deduce the estimates 



a = 0.847(3) and (f> = 0.867(3) 



(2.11) 



Importantly this is consistent with the estimates found in 
|13j for the growth process on the square lattice where 4> 
was estimated to be 0.88(7), subsequently confirmed in- 
dependently and directly for the equilibrium model in 
[T5] where — 0.84(3) was found. Our estimate of 
acf) = 0.734(6) implies that the exponent involved in 
the scaling of the ej, namely (1 — «)</>, takes the value 
0.133(4). In Figure [4] we plot 62(71) and 63(71) against 
We immediately see the strong corrections to 
the scaling forms still apparent even at the long lengths 
we have simulated. The maximum in 62(71) at 71 ss 1000 
mirrors the behaviour found in /2 (71) in Figure [3j 

In Figure |4] we have also marked the asymptotic val- 
ues E2 and -B3 which can be deduced from the following 
argument. Consider the case when a trail has formed a 
large n-step loop which occupies 771 lattice sites. Any site 
of this loop could have been the starting point. In order 
for this site to be visited only once, the loop must have 



closed at the first return visit, which occurs with a proba- 
bility of (|). In order for this site to be visited twice, the 
loop must have closed at the second return visit, which 
occurs with a probability (1— |)(|). Finally, for this site 
to be visited three times, the loop closes at the third re- 
turn visit, which occurs with a probability (1 — |)(1 — |). 
Therefore we find for large loops the asymptotic values 
mi/m = 1/5, 7712/771 = 4/15, and m^/m = 8/15. Using 
(2.1 ), we obtain 



61(71) 



62(71) 



35' — 35' ^^^")^35 

(2.12) 

as 77 — >■ c», whence we identify E2 = 4/35 « 0.114 and 
E3 = 8/35 w 0.228. 
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FIG. 4: Plot of the average numbers of doubly, 62(71), and 
triply, 63(11), visited sites against l/n"'^"^^. Their asymptotic 
values, E2 = 4/35 ^ 0.114 and E3 = 8/35 « 0.228 respec- 
tively are marked as filled circles. 



III. CANONICAL ISAT MODEL ON THE 
TRIANGULAR LATTICE 

The canonical model [12] of self-interacting trails 
(ISAT) on the triangular lattice is defined as follows. 
Consider all different bond-avoiding paths (open trails 
and loops) ipn of length 71 that can be formed on the 
triangular lattice with one end fixed at a particular site 
(the set ri„). Associate an energy — e with each doubly- 
visited site and an energy — 2e with each triply-visited 
site. For each configuration ipn count the number m2{^Pn) 
of doubly- visited sites and 7713 ((/3„) of triply- visited sites 
of the lattice and give that configuration a Boltzmann 
weight fjj'^^+'irrii ^ where lu = exp (/?£). The partition func- 
tion of the ISAT model is then given by 



E 



UJ 



(3.1) 



We use the superscript (2) to denote the canonical model. 
The choice of this notation will become clear in the next 
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section. The average of any quantity Q over the ensemble 
set fi„ of allowed paths of length n is given generically 

by 



(Q)n(^) 



(3.2) 



The reduced free energy {uj) per step is given by 

4')M = ilogZ(2)(c.), (3.3) 
n 

and the internal energy u^\uj) and specific heat c"n\i^) 
can be found as the first and second derivatives of the 
free energy with respect to /3. We will also require the 
third derivative tn (w) . However, these can all be found 
from the expected values of moments of m2 and m^. The 
internal energy (setting e = 1) is given by 



,(2) 



{1712 + 2to3 



and the (reduced) specific heat is given by 



{{1712 + 2to3)^) ~ ((^2 + 2m3))^ 



(3.4) 



(3.5) 



We have simulated the canonical ISAT model using 
FlatPERM [24] for lengths up to 1024, generating 1.6 x 
10*" samples at that length. First we consider a high 
temperature point at w = 1.5 and show that the scaling 
of the size of the trail is in accord with = 3/4 and 
the standard swollen phase expectations: in Figure [S] we 
display a plot of log i?^ against log n. 



of 




FIG. 5: A plot of log_R„ against logn at {u}2,uj3) = (1.5,2.25) 
in the swollen phase showing a slope of 2v — 3/2. 

By considering the specific heat we find a weak phase 
transition in contrast to that found on the square lattice. 
There is little sign that the specific heat diverges: in 
Figure [6] the value of the maximum of the specific heat 
is plotted against log n. 




FIG. 6: Plot of the value of the maximum of the specific heat 

(2) 

gainst logn. This suggests that the specific 
heat does not diverge as the polymer length is increased. 



We recall that the 6'-point transition of ISAW has this 
non-divergent specific heat behaviour, and so we have 
considered the third derivative of the free energy with 

(2) 

respect to /3, t„ (w), which would diverge slowly for the 
6'-point. The third derivative of the free energy is the first 
derivative of the specific heat. In Figure [7| we display the 
absolute value of the two peaks of the third free energy 
derivative. The third derivative has two peaks: one posi- 
tive and negative in value. They show a weak divergence: 
we have extracted local exponents in both cases, which 
would be given by (1 -I- a)4>, were this a critical point. 
We find the values 0.23(6) and 0.35(6): this is consistent 
with the ISAW fl-point, where (1 + a)(f> = 2/7 « 0.28. 
While there are clearly very strong corrections to scaling 
still in the data, no divergence is found in the specific 
heat and a weak one in the third derivative. Therefore it 
is tempting to conjecture that the canonical ISAT model 
on the triangular lattice has a collapse transition that 
lies in the 0-point universality class, rather than square 
lattice ISAT collapse universality class. 



IV. EXTENDED ISAT MODEL ON THE 
TRIANGULAR LATTICE 

A. The model 

The extended model of self-interacting trails (elSAT) 
on the triangular lattice is defined as follows. Consider 
all different bond-avoiding paths ipn (open trails and 
loops) of length n that can be formed on the triangu- 
lar lattice with one end fixed at a particular site (the set 
r2„). Associate an energy —£2 with each doubly- visited 
site and a different energy —£3 with each triply-visited 
site. For each configuration ipn count the number m2{fn) 
of doubly- visited sites and m3{ipn) of triply-visited sites 
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FIG. 7: Plot of the height of the peaks of t^n^uj), the 
third derivative of the free energy with respect to temper- 
ature against n. The third derivative has two peaks: one 
positive and one negative in value. The top figure shows 

(2) (2) 

t„ — maxuj tn and the bottom figure shows t„ — min„ t„ . 



of the lattice and give that configuration a Boltzmann 
weight Lu^P^Lu^^'', where ujj — exp{f3£j). The partition 
function of the elSAT model is then given by 



(4.1) 



The probability of a configuration (p„ in the equilibrium 
model is 



Pe 



(^„;c.2,W3) = ojT"^''"^^r^^"^/ZnioJ2,oj,) . (4.2) 



The average of any quantity Q over the ensemble set of 
allowed paths fin of length n is given generically by 

J2 Qi^n) PEiiPn) ■ (4.3) 

Let us define the expectation values of the number of j- 
fold visited sites per unit length Uj for the equilibrium 



model as 



tf>(n) 



and the fluctuations in these as c,- with 



c-,(n; W2, W3) 



(4.4) 



(4.5) 



Note that from Equation (2.1 1, which holds for trail con 



figuration regardless of how they are generated, we have 
ui + 2u2 + 3u3 ^ 1 . (4.6) 
The canonical ISAT model is then given by the restriction 



UJ2 



and 



UJ3 



(4.7) 



By fixing the ratio of the energies 63/82 = k we have the 
generalisation 



W2 



and 



UJ3 



(4.8) 



which gives a family parametrised by k of generalised 
one-parameter ISAT models. This gives the partition 
function of our fixed k versions of the elSAT model as 



Setting £2 = 1 the internal energy is 



^W(^)^ ((^2 + fcm3)) 



and the (reduced) specific heat is 



((m2 + kms)"^) - ((to2 + km^))"^ 



(4.9) 



(4.10) 



(4.11) 



The canonical ISAT model has fc = 2. 

We have simulated the general elSAT with a two- 
parameter FlatPERM algorithm [21] up to lengths 128, 
generating 8.8.6 x 10^ samples at that length. We have 
also simulated various specific sub-cases for fixed values 
of k via one-parameter FlatPERM algorithm typically 
up to length 1024, generating in each case roughly 10^ 
samples at that length. 



B. Mapping of the growth process to the 
equilibrium model 

The growth process on the square lattice has been 
mapped [52 [23] to a specific Boltzmann weight of the 
equilibrium model by considering loops. The probability 
distribution Pcifn) of the growth process on the trian- 
gular lattice can be rewritten in terms of mj{ipn) as 



Pci^n) = 



nil (Vn)-l 



ri-1 



1 

15 



m2 {^n ) 



1 

15 



(4.12) 
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using the relationship given in Equation (2.1 1. We then 
notice that by setting uii = 5/3 and ^3 = 25/3 the equi- 
Ubrium model has probability distribution 



Ve 



5 25 



1 



3' 3 

and so we can deduce 



7 (i £5^ 
"^"^3' 3 



m2(ipn) 



5 25 
3' "3" 



(4.13) 



(4.14) 



Note that the normalisation is different though since the 
sum over all walks of fixed length n gives the probability 
of walks being still open in the case of the growth process, 
and unity in the case of the equilibrium model. 

We now consider the probability of producing configu- 
rations in the growth process conditioned on the process 
having continued to length n, that is, we consider 



which implies that 



Hence we have 



PE \(P' 



5 25 

3' "3" 



noting that this implies that 

Z„(^,|^=65"-ip„ 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



In any case, when we simulate the growth process we are 
effectively simulating the equilibrium elSAT model at the 
point (0^2,1^3) — (5/3,25/3). We immediately have that 



ej{n) 



5 25\ , , . ^ / 5 25 

n; - , — I and fAn) ^ c, [ n 



3' 3 
(4.19) 

Using the definition of the one-parameter family of mod- 
els via Equation (4.81, the growth process is equivalent 
to w = 5/3 in the model where 

log(5/3) 

Importantly, this is not the value of k for the canonical 
ISAT model (with fc = 2) on the triangular lattice that 
has been studied via exact enumeration 1121. 



C. The k = kc elSAT model 

As described in Section 2 we have simulated the growth 
process and found divergent fluctuations in the number 



of doubly and triply-visited sites — a sign of critical be- 
haviour. Now that we have mapped the growth process 
onto a specific temperature of the k = kc elSAT model 
we can verify that this point is indeed the collapse tran- 
sition point. We have simulated the k = kc model up to 
length 1024 using a one-parameter FlatPERM algorithm. 
We find a divergent specific heat with a single pronounced 
peak near uj = 5/3. We begin by finding the location and 
size of the peak of the specific heat of the model. In Fig- 
ure [8] we plot the logarithm of the peak height of the spe- 
cific heat against log n along with a line corresponding to 
the exponent value ai/) — 0.734, which was obtained from 
the growth process. Moreover, using the finite length es- 
timate of (f> — 0.867 obtained from the growth process to 
extrapolate the location of the transition (see Figure [9]) 
we find an estimate of ojt = 1.669(4): this is consistent 
with the growth process point at uj — 1.66 being the 
transition point. 
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FIG. 8: Plot of the logarithm of c„ = maxi^ <^n^^ ■, the value of 
the maximum of the specific heat, against logn. The straight 
line has slope 0.734. 



While the exponent estimates clearly discount a first- 
order transition, given that there was a question about 
the first-order nature of the ISAT model on the square 
lattice at the kinetic growth point we now show how the 
distribution of the numbers m = -f kcms changes as 
the temperature is moved through the transition point: 
this is displayed in Figure [TO] There is no sign of first- 
order behaviour in the distributions. 

To complete our numerical analysis we provide a scal- 
ing plot around the transition of the specific heat using 
the exponents values from the growth process. This is 



found in Figure 11 and it shows an excellent fit to the 
crossover scaling form |19] 



n°"^C{{uj - UJt)n''') . 



(4.21) 



Our data at length 1024 there fully confirms that the 
kinetic growth point is the location of a phase transition 
and, moreover, that this transition is critical in nature. 
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FIG. 9: Plot of the location, ujmax, of the peak of the specific 
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FIG. 11: Scaling plot of the specific heat around the transition 
temperature, using the exponents from the growth process. 



D. The fc = elSAT model 




To explore this fc-dependency more we now examine 
the A; = elSAT model. When fc = only doubly-visited 
sites are given a Boltzmann weight so that 



The internal energy is 



,(0) 



and the specific heat is 



.(0) 



(^) 



{m2) 



(4.22) 



(4.23) 



(4.24) 



FIG. 10; Plot of the distribution p„(m/n), where m = 
m2 + kcm^, at temperatures near, and at, the tempera- 
ture at which the specific heat attains its maximum for 
length n — 1024. The specific heat attains its maximum at 
uj — uJmax = 1-69 and the distribution is plotted for this value 
and at oj = 1.63, 1.66, 1.72, 1.76: the plots move from left to 
right as oj is increased. 



The exponents values indicate that the ISAT model on 
the square lattice and the k = kc elSAT model on the 
triangular lattice are in the same universality class. This 
is in contrast to the observation above that the canoni- 
cal ISAT models on the two lattices are not in the same 
universality class. This inexorably leads us to the con- 
clusion that the universality class of the phase transition 
of the elSAT model on the triangular lattice depends on 
the value of k. 



In Figure [12] the value of the maximum of the specific 
heat over all ui is plotted against logn. In stark contrast 
to the behaviour of the k = kc model but in common 
with the k — 2 canonical model the specific heat does 
not seem to diverge as the length is increased. In com- 
mon with the approach we took for the k = 2 model, we 
have also examined the third derivative of the free energy. 
In Figure [13] we see that the absolute maximum of this 
quantity is weakly divergent with an exponent of 0.23(6). 
This compares well to the 6'-point value of 2/7 « 0.28, 
especially given the relatively short lengths of the sim- 
ulations. It is tempting then to conjecture that fc = 
and k — 2 models have collapse transitions in the same 
universality class and that this class is the 0-point. 



E. The Triple model ("fc = oo") 

Next we consider a model where only triply-visited 
sites are weighted — essentially a "fc — oo" elSAT — 



9 




Note that only in the Triple model u) 
finite k models we denote uj = uj2- 



UJ3. 



while for the 



As a function of w we find evidence for a very strong 
phase transition. In fact when we considered the scaling 
of the peaks of the specific heat they scale faster than 
linearly (linearly is the maximum theoretical asymptotic 
behaviour). Since linear behaviour would indicate a first- 
order phase transition we considered the distribution of 
triply-visited sites at temperatures around that which 
gives rise to the peak of the specific heat: this is plot- 
ted in Figure 14 The distribution displays an unam- 



0.06 



FIG. 12: Plot of the value of the maximum of the specific heat 
gainst logn. This suggests that the specific 
heat does not diverge as the polymer length is increased, as 
is the case in the canonical model (fc = 2). 
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FIG. 13: Plot of the height of one of the peaks of tli\uj), the 
third derivative of the free energy with respect to temperature 
against n. The third derivative has two peaks: one positive 
and one negative in value. The figure shows tn = min^^ ti"' 
which is the larger in absolute value and the one seemingly 
less affected by corrections to scaling. 



that is, the partition function is given by 
The internal energy is 

n 

and the specific heat is 

" n 



(4.25) 



(4.26) 
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FIG. 14: Plot of the distribution pn(m^/n) of triply- visited 
sites for the Triple model at temperatures near, and at, the 
temperature at which the specific heat attains its maximum 
for length n = 1024. The specific heat attains its maximum 
at = ujrnax = 7.41 aud the distribution is plotted for this 
value and a,t u) = 7.31,7.34,7.48,7.52: the plots move from 
left to right as cj is increased. 



biguous double peaked form which is a clear sign of a 
first-order phase transition. Note the sharpness of the 
transition in temperature and how the distribution moves 
from being peaked around low values of triply-visited 
sites just below the transition temperature to peaked 
around large values just above the transition tempera- 
ture. This is classic first-order behaviour. To get an 
idea of where the transition takes place in the thermody- 
namic limit we have extrapolated the transition location 
in Figure |15| We estimate a thermodynamic limit value 
of LOt{Triple) = 6.96(6). 

We now see that when triply-visited sites are weighted 
the transition is clearly first-order, while when doubly- 
visited sites are the only ones given additional weight 
the transition is a weak second-order transition, not un- 
like the 6'-point. In between, the growth model displays 
intermediate behaviour. 



10 




0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 0.004 
1/n 



FIG. 15: Plot of the location, ujmax, of the peak of the specific 
heat against 1/n for the Triple model. 



F. General elSAT model 

To understand whether the first-order nature of the 
collapse transition persists when doubly-weighted sites 
are given some weight we have examined the fc = 6 model. 
Once again the divergence of the specific heat is very 
strong, and plotting the distribution of triply- visited (see 
Figure 16) the classic doubly peaked form is once again 
apparent. It seems that the first-order nature persists 



0.0035 




FIG. 16: Plot of the distribution p„(m3/n) of triply- visited 
sites for the k = 6 model at temperatures near, and at, the 
temperature at which the specific heat attains its maximum 
for length n = 1024. The specific heat attains its maximum 
a.t uj — uJmax = 1.419 and the distribution is plotted for this 
value and at a; = 1.410, 1.414, 1.425, 1.430: the plots move 
from left to right as u is increased. 

for large finite values of k. Given the evidence available, 
the simplest scenario that presents itself is the following: 
For values of k > kc the elSAT model on the triangular 
lattice displays a first-order phase transition, while for 



k < ka the model displays a weak second-order transition 
presumably in the universality class of the IS AW 0-point. 
Separating these two behaviours is the k = kc model, 
where the transition is second-order but with a divergent 
specific heat in the universality class of square lattice 
ISAT. This implies that the growth process point is a 
multi-critical point that is the meeting of a line of first- 
order transitions to a line of critical ones. 

The question then arises as to the nature of the low 
temperature phase for different k. In Figure 17 we give 



a density plot of the largest eigenvalue of the matrix of 
second derivative of the free energy with respect to the 
two variable uj2 and W3 at length n = 128: this allows 
us to search for any other possible transitions. Intrigu- 
ingly this plot show evidence for transitions at large uj2 
and ^3. This would seem to indicate a collapse-collapse 
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FIG. 17: Density plot of the logarithm of the largest eigen- 
value Xrnax of the matrix of second derivatives of the free 
energy with respect to CJ2 and UJ3 at length n = 128 (the 
lighter the shade, the larger the value). 

phase transition with two different collapsed states. Such 
a transition has been seen to occur in the low temperature 
behaviour of the ISAW when a stiffness energy is added 
to model semi-flexible polymers. One of the two low tem- 
perature phases in that model is an ordered crystalline- 
like phase. Relatedly this occurs when the collapse is a 
first-order transition. 

To test the possibility that the low temperature col- 
lapsed phase when k > kc is frozen we have plotted in 
Figure [18] the proportion of steps of the trail not involved 
with a triply-visited site. Given the collapsed nature of 
the phase it is appropriate to assume that the polymer 
has a surface [25] and so we plot our estimate against 
1/n^/^. We have chosen a point in the fc = 6 model. The 
quantity [1 ~ 3^3(71)] is tending to zero as n is increased, 
suggesting that, proportionally, all sites of the trail will 
be triply-visited: we can deduce that configurations pro- 
duced are maximally dense and solid-like. In Figure [19] 
we display a typical configuration produced by our simu- 
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FIG. 18: Plot of 1 — 3u3{n), which measures the proportion 
of steps that are not involved with triply- visited sites per unit 
length, against l/\/n at a point (1.58, 15.6) in the hypothe- 
sised frozen (crystal-like) phase. As the length increases this 
quantity vanishes. 



FIG. 20: Plot of 1 — 3M3(n), which measures the proportion 
of steps that are not involved with triply- visited sites per unit 
length, against l/\/n at a point {uj2,uj3) = (4, 16) in the col- 
lapsed liquid-drop-like globule phase. As the length increases 
this reaches a non-zero value. 




FIG. 19: A typical configuration at length 512 produced at 
(aj2,i^3) = (1, 10) which looks like an ordered crystal. 



lations when (a;2, W3) = (1, 10): it is space-filling, forming 
a crystal-like structure, being almost totally made up of 
triply-visited sites. We find the same behaviour for the 
k ~ ka model at low temperatures. 

For the sake of comparison we have considered the 
same quantity 1 — 8^3 (n) in the low temperature region 
of the k = 2 model. In Figure [20] we plot the proportion 
of steps that are not involved with triply-visited sites per 
unit length against l/n^/^ at the point (a;2, W3) = (4, 16): 
this quantity clearly converges to a non-zero value. In 



Figure 21 we display a typical configuration produced 
by our simulations when (w2,'^3) — (5,1): while dense, 
the configuration still has internal holes and seems disor- 
dered. 

From the above considerations it seems that the phases 



FIG. 21: A typical configuration at length 512 produced 
at (102,103) = (5,1), which is in the globule phase: it looks 
disordered and rather more like a liquid-like globule than a 
crystal. 



that exist in our elSAT model on the triangular lat- 
tice are of a similar type to those in the semi-flexible 
ISAW phase diagram [5HTU] : namely a high temperature 
swollen phase and two low temperatures phases: one a 
disordered liquid drop-like globular phase and a second 
crystalline-like phase. 



V. CONCLUSIONS 

We can now build a complete picture of the phase dia- 
gram of the extended ISAT model on the triangular lat- 
tice. A schematic of the conjectured phase diagram is 
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shown in Figure [22] For small U2 and ^3 we see the usual 
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FIG. 22: Schematic of the proposed phase diagram of the 
extended ISAT model on the triangular lattice. The filled 
circle is at the location of the kinetic growth point, and the 
open circles represent estimates of the collapse transition for 
various values of k. 

swollen polymer phase where v = 3/4. For large enough 
L02 regardless of wa we find a collapse phase as occurs 
in the ISAW model and a transition between the swollen 
and collapsed globule phases which seems to be 0-like. On 
the other hand for large enough u)^ we find crystal-like 
configurations that are space filling and internally con- 
tain only triply-visited sites. Between the swollen phase 
and the crystal-like phase the collapse transition is first- 
order. Separating this line of first-order transitions from 
the line of 0-like transitions is a multi-critical point. This 
point is precisely the point to which the kinetic growth 
process of trails maps. 



It would seem that this larger parameter space has ex- 
posed a way of understanding the apparent differences 
and similarities of the ISAT and ISAW models. When 
ISAW is generalised with the addition of stiffness, and 
ISAT on a large enough coordination number lattice is 
generalised with different weightings for different num- 
bers of visits, both display three polymer phases: swollen 
coil, globule and polymer crystal. Note that for the 
elSAT model considered here, the crystal-like phase has 
non-zero entropy and as such is strictly speaking not a 
proper crystalline phase, but rather a maximally dense 
phase. This is in contrast to the polymer crystal phase 
found in the semi-flexible ISAW model [5HTU] , where one 
finds zero-entropy /?-shcet like structures. 

In two dimensions both semi-flexible ISAW and elSAT 
have a first-order collapse transition between swollen coil 
and crystal-like phases, and a 0-point like transition be- 
tween swollen coil and globule phases. It is also apparent 
that the square lattice ISAT model is unusual in that it 
only displays the multi-critical point, which is found in 
these generalised models as part of a larger phase dia- 
gram. 
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